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ABSTRACT 

We report observations of 47 candidate protostars in two optically thick lines 
(H2CO 2i2 - 111 and CS 2 - 1) and one optically thin line (N2H+ 1-0) using the 
IRAM 30-m, SEST 15-m, and Haystack 37-m radio telescopes. The sources were 
selected from the redness of their spectra (Tboi < 200 K) and their near distance 
{d < 400 pc). Most of the sources have asymmetric optically thick lines. The observed 
distribution of velocity differences, dV = {Vthick — ythin) / ^ythim is skewed toward 
negative (blue-shifted) velocities for both the H2CO and CS samples. This excess 
is much more significant for Class than for Class I sources, suggesting that we 
detect infall motions toward Class and not toward Class I sources. This indicates a 
difference in the physical conditions in the circumstellar envelopes around Class I and 
Class sources, but does not rule out the presence of infall onto Class I sources by e.g. 
lower opacity gas. Bipolar outflows alone, or rotation alone, cannot reproduce these 
statistics if the sample of sources has randomly oriented symmetry axes. We identify 
15 spectroscopic infall candidates, of which 6 are new. Most of these infall candidates 
have primarily turbulent rather than thermal motions, and are associated with clusters 
rather than being isolated. 



Subject headings: Stars:formation, ISM:kinematics and dynamics, ISM:molecules 



1. Introduction 

Low mass stars form by gravitational collapse of dense cores in molecular clouds. The 
collapse is believed to proceed from a condensed initial state (Larson 1969, Shu 1977) with a rapid 
development of a system consisting of a protostellar core which accretes mass from a disk and 
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envelope. This gravitational collapse model can now be explored through detailed observations 
and numerical calculations. Observational studies of star-forming infall are indispensable to 
understand the kinematics of gravitational collapse onto young stellar objects (YSOs). This, 
in turn, is essential to study the origin of the initial mass function and of multiple systems. 
Observations are now beginning to probe a variety of size scales, density and temperature regimes 
in many sources. 

In recent years, millimeter- wavelength observations have increased the evidence for infall 
motions onto nearby YSOs. The main procedure is to observe spectral lines that trace high 
densities (n > 10^ cm~^) and have red-shifted self-absorption spatially concentrated around an 
embedded YSO. The spatial concentration and the use of dense gas tracers are important to 
ensure that the observed kinematic signature is associated with the YSO. For instance, Zhou et 
al. (1993) mapped the globule B335 in the 2i2 — In and 3i2 — 2ii optically thick lines of H2CO 
and in the optically thin lines of C^^O 1 — and C^^S 2 — 1. They reproduced the line profiles on 
spatial scales of 0.02 pc from the source with an inside-out collapse model (Shu 1977). Choi et al. 
(1995) modeled the same B335 observations with a Monte Carlo radiative transfer code (Bernes 
1979), confirming the numerical results from Zhou et al. (1993). Later Zhou (1995) modeled the 
observed H2CO and CS line profiles toward IRAS 16293-2422 as arising from a combination of 
infall and rotation, using the Terebey, Shu & Cassen (1984) collapse model. The case for infall 
motions in these sources is strong because of the detailed observations and modeling that has been 
done on them. However, in such a small sample of candidates, peculiarities of source structure 
and kinematics (eg, outflows or rotation) could possibly account for the apparent signatures of 
infall. Therefore, these well studied cases are not sufficient by themselves to support a claim for 
observed infall motions onto a wide class of sources. 

Another approach to search for infall evidence is to observe a statisticalh^ large sample 
of sources in two lines tracing high density gas, one optically thick and the other optically 
thin. In such a sample, rotation and bipolar outflow motions having symmetry axes in random 
orientations would tend to produce optically thick line profiles with equal numbers of red-shifted 
and blue-shifted self-absorption. On the other hand, infall motions in centrally concentrated 
regions would only produce red-shifted self-absorption in the line profiles (Leung &; Brown 1977). 
A spectral line survey of a large source sample should reveal whether there is a statistically 
significant excess of sources with red-shifted or blue-shifted self-absorption. This approach docs 
not depend much on the details of the models, but mainly on whether there is a prevalence of 
inward or outward motions. At the same time such a study can yield a list of collapse candidates 
for further detailed studies. Recent surveys aimed to find candidate collapsing YSOs include 
Mardones et al. (1994) which used the C3H2 2i2 — loi line; Wang et al. (1995) who observed the 
H2CO 3i2 — 2ii line in a sample of Bok globules, and Gregersen et al. (1997) which used HCO^ 
4-3 and 3-2 lines in a sample of Class sources (Andre, Ward-Thompson Sz Barsony 1993). None 
of these surveys has emphasized the statistical properties of the samples because of the limited 
number of targets observed. 
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In this paper we present the results of a molecular line survey toward 47 of the reddest nearby 
YSOs, more than twice the number of sources in any of the previous surveys. In section 2 we 
discuss our source selection criteria, the molecular lines and the telescopes we used. In section 3 
we describe the normalized velocity difference (SV) between the peaks of optically thick and thin 
lines, as a quantitative measure of the line profiles. We then analyze statistically the distribution 
of SV for different molecular lines and source subsamples. In section 4 we interpret the statistics 
kinematically. Finally, we summarize our conclusions in section 5. 

2. Observations 

2.1. Source Sample 

The sources were selected solely based on their observed spectral energy distribution (SED) 
and estimated distance from the Sun in order to be unbiased with respect to prior knowledge 
of their molecular line shapes. To select the youngest sources (best infall candidates) we used 
the bolometric temperature (T{,o^; Myers, &i Ladd 1993) whenever known, otherwise estimated 
T{)oi based on IRAS and submm continuum observations. T{,oi is defined as the temperature 
of a black body with the same mean frequency as the source spectrum. Chen et al. (1995) 
found that Class sources have T^oi < 70 K, and Class I sources (Lada &: Wilking 1984) have 
70 < Tjjoi < 650 K. We imposed a limit of Tj,gi < 200 K to select the most embedded Class I and 
Class sources. We relied mostly on the papers by Ladd, Lada, & Myers (1993) in Perseus and 
Chen et al. (1995, 1996) in the Taurus, Ophiuchus, Chamaeleon, Lupus and Corona Australis star 
forming clouds for determination of T{,oi. Chen et al. required the source fluxes to be known in at 
least 6 different wavelengths in order to compute accurate values of Tfj^i. That condition is easily 
met by most IRAS sources which have also been detected at optical or near-infrared wavelengths. 
However, many of the most embedded Class I sources and most Class sources do not satisfy that 
requirement. Only a handful of Class sources have been observed in enough different mm and 
submm wavelengths to make a good determination of Thd possible (eg. B335 and L1527, Ladd et 
al. 1991). We therefore added sources to our list based on more relaxed criteria in order to have a 
larger sample of the youngest known nearby YSOs. 

We required IRAS point sources in our sample to have rising spectra between 25 and 100 fim 
(with flux ratios 2-4 for wavelength pairs of 25/12, 60/25 and 100/60/7,m), and to be located close 
to the centers of their parent molecular cloud core maps. These conditions are necessary to exclude 
galactic cirrus and extragalactic objects which may have high 100/60)um flux ratios. We also 
added sources with known bright mm or submm continuum emission, and nearby Class sources 
known to us by mid 1995. We consulted Casali et al. (1992) for submm continuum observations in 
Serpens, and Henning et al. (1993) and Rcipurth et al. (1993) for submm continuum observations 
toward southern PMS stars. We also consulted Persi et al. (1990), Carvallo et al. (1992) and 
Bourke et al. (1995) for IRAS and near-infrared observations of southern PMS stars, many of 
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them located in Bok globules. We determined T^oi for all the sources detected in at least three 
different wavelengths to include sources not quoted in Chen et al. (1995, 1996). The derived values 
of Thoi for sources with poor wavelength coverage are upper limits if they lack long wavelength 
fluxes, and lower limits if they lack short wavelength fluxes. Because of our selection criteria, 
we are confident that Thoi < 200 K in all sources in our sample. To determine a representative 
systematic error in T^o; in sources with poor short wavelength spectral coverage, we determined 
Tboi including and excluding wavelengths shorter than 60 /xm in sources with good spectral 
coverage. We found that the error in Tf^oi when neglecting wavelengths shorter than 60 fim is 
typically < 20 K for sources having T^oi < 100 K. Our determination of T^oi is consistent with 
the values given by Gregersen et al. (1997) also with limited wavelength coverage. 

Table || lists all the sources detected in at least one optically thin and one optically thick line. 
Columns 1 and 2 identify the source, columns 3 and 4 list the observed coordinates (B1950.0), 
column 5 lists the distance from the Sun in pc, column 6 lists Tf,oi, and column 7 the number 
of continuum fluxes used to derive T^o;. We give lower or upper limits to T^oi rather than values 
for the sources IRAS 03282+3035, HH211mm, VLA1623, and S68N due to lack of wavelength 
coverage close to the peak of their spectral energy distributions. Although these estimates are 
crude, they are sufficient for our analysis in section 3, since we only need to know whether T^oi is 
below 70 K or not. Column 8 lists the references we consulted to determine values of Tboi- Finally, 
column 9 indicates whether there is outflow evidence from the source either as high velocity 
molecular line wings (CO) or as shock fronts (HH). 

Even though we did not use the outflow properties as a source selection criterion, 87% of the 
sources in Table || actually show evidence of driving an outflow. This is expected when selecting 
the most embedded sources, since they tend to drive strong molecular outflows (Bontemps et al. 
1996). YSO outflow models usually require mass infall onto an accretion disk (see references in 
Bachiller 1996), so there is reason to expect that most sources in this survey must have infalling 
envelopes. However, we emphasize that in this paper we are seeking to examine direct kinematic 
evidence for infall motions onto the sources. 

We did not detect 19 of the sources from our initial target list in either the H2CO line, the 
N2H+ line, or both. Table || lists the positions and rms noise (main beam brightness temperature) 
of the undetected lines for reference. We detected the N2H"'" line toward the sources IRAS 
11436-6017 and GSS30-1, but it was too weak for a good velocity and width determination (peak 
Tmb of 0.15 and 0.3 K respectively). We note that the average T},oi of these sources in Table |2| is 
higher than that of the detected sources listed in Table |l|. In fact, sources with T^oi < 50 K in 
Table | were easily detected in all the observed lines. 
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2.2. Lines 

It is important to observe both optically thick and optically thin lines in order to distinguish 
double-peak profiles produced by optical depth effects from those due to two velocity components 
along the line of sight. The observed molecules should have a high dipole moment in order to be 
sensitive to the high density gas in the immediate neighborhood of the YSO and not to the more 
extended, lower density gas surrounding the parent cores. Also, the molecular species should be 
abundant in molecular cloud cores, but unaffected by peculiar chemistry effects. For example, SiO 
and CH3OH are well known for their enhanced abundances by many orders of magnitude in the 
outflow lobes of some YSOs (e.g. Bachiller 1996), presumably due to desorption of molecules from 
grain mantles in shock fronts. 

Following these criteria, we chose the CS {fi = 1.96 Debye) and H2CO (/x = 2.33 Debye) 
molecules to carry out the survey observations. H2CO was our first choice molecule, following the 
example of Zhou et al. (1993). The abundance of H2CO is known to be only moderately enhanced 
in bipolar outflows (Bachiller, & Perez 1997). We subsequently observed part of the survey in the 
CS molecule in order to check that our results were not being biased by peculiar chemistry in 
low mass star forming regions. We believe the CS and H2CO lines are generally optically thick 
because of the observed brightness temperatures and line shapes. In a few cases, when the line 
is weak {Tmb ^ 1 K) and approximately of Gaussian shape, it is likely to be optically thin (eg. 
L673A, 04325+2402, L1172; see figs. II). 

Ideally, one would observe an isotope of the optically thick tracer molecule as the optically 
thin tracer. Then both lines would probe the same density, temperature and chemistry, and 
would be observed with essentially the same angular resolution. However, the 2-mm H2^C0 and 
the 3-mm C^'^S lines are in general too weak for accurate determination of line velocities in a 
large source sample. We chose to observe the N2H"'" J=l — line instead (fi = 3.4 Debye). This 
molecule has seven hyperfine components, in two groups of three lines and one isolated component 
{JF1F2 = 101 — 012). Caselli, Myers, & Thaddeus (1995) determined the frequencies of all 
components with a precision of 6 kHz (0.02 km s~^). The hyperfine structure of N2H^ also permits 
a determination of the optical depth in the line, thus, allowing an internal check on the thinness 
of the isolated N2H"'" 101 — 012 line. Caselli et al. showed that among the seven components, only 
two show any evidence of anomalous excitation {F1F2 = 12 — 12 and F1F2 = 10 — 11), and that 
these anomalies represent a departure of < 20% from the value expected in LTE. Any such small 
anomaly in hyperfine line ratios is included in the estimated uncertainty of the optical depth. 
This uncertainty is generally much smaller than the estimated optical depths. Womack, Ziurys, 
& Wyckoff (1992) found that the N2H"'" J= 1 — line is narrow in dark clouds, with no evidence 
of line wings. The N2H"*" line- widths and hyperfine components make it an ideal probe of the 
systemic velocity toward young stellar objects. 
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2.3. Telescopes 

We used the IRAM 30-m telescope in Pico Veleta, Spain, in 1995, May 10-12. We observed 
simultaneously in the 2 and 3mm atmospheric windows with receivers tuned to the ortho H2CO 
2i2 - 111 (140.839518 GHz) and N2H+ 101 - 012 (93.176265 GHz) lines respectively. The line rest 
frequencies were taken from Lovas (1991) and Caselli et al (1995) respectively. The backend had a 
digital correlator with spectral resolutions of 20 and 10 kHz at 2 and 3mm wavelengths respectively 
(0.04 and 0.03 km s~^). We used overlap frequency switching by ±3.85 MHz because this method 
is efficient and ensures that the observed line profiles are free of emission from an off-position. We 
checked pointing regularly with continuum sources and found it good to better than 5". Typical 
system temperatures were ~400 K at 2-mm and ^^350 K at 3-mm. The main-beam efficiency at 
93 and 141 GHz was taken to be 73% and 55% respectively. 

We used the SEST 15-m telescope in La Silla, Chile, in 1995, September 19-20 and September 
24, to observe simultaneously the same 2-mm and 3-mm H2CO and N2H^ lines as with the IRAM 
30-m telescope. The weather was excellent, with system temperatures around 150 K in both lines. 
We splitted the high resolution acousto-optical spectrometer into two bands, providing a spectral 
resolution of 43 kHz in each band (0.138 km s'^at 93 GHz and 0.11 km s'^at 141 GHz). The 
acousto-optical spectrometer sometimes produces spurious emission features in the central 3-5 
channels; so to avoid interference with the observed spectral lines, we centered the spectrometer 
off the rest velocity of each source by 2-8 km s~^. We also observed a subsample of sources in the 
para H2CO 2o2 - In transition at 145.602953 GHz and the C^^S 2-1 transition at 96.412982 
GHz. We used main beam efficiencies of 74% at 93 GHz and 96 GHz, 68% at 141 GHz, and 67% 
at 145 GHz as indicated in the SEST handbook. 

In 1995, December 7-11 and 29-31, we used the NEROC Haystack]^ 37-m telescope to observe 
the sources visible from Haystack in the CS 2-1 line at 97.980968 GHz. We observed left and right 
circular polarizations simultaneously with two SIS receivers. The backend had a digital correlator 
which we used with 17.8 MHz bandwidth and 13 kHz spectral resolution (0.04 km s^^). Typical 
system temperatures were 200 K. The antenna main beam FWHM at 98 GHz (21") compares well 
with the IRAM 2-mm and 3-mm beams. We measured a main beam efficiency of ~ 16% it 3% by 
scaling spectra of IRC+10216 to those taken at the IRAM 30-m telescope by Mauersberger et al. 
(1989). 

3. Results 



^Radio astronomy observations at the NEROC Haystack Observatory of the Northeast Radio Observatory 
Corporation are supported by a grant from the National Science Foundation. 
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3.1. Variety of Observed Spectral Line Profiles 

Figure |^ shows the observed spectra toward all the Class sources {Tboi < 70 K) in our 
sample ordered by right ascension, and figure |2| the spectra obtained toward the Class I sources 
(70 < Tboi < 200 K). The N2H+ 101 - 012 line is drawn at the bottom of each box; H2CO 
2i2 — 111 in the middle, and CS 2 — 1 at the top if it was observed. In a few sources, such as IRAS 
16244-2422, the N2H"*" 101 — 012 line is very weak, but the velocity (dashed line) is sufficiently well 
determined because it is based on all seven hyperfine components, not just the one shown. We find 
a wide variety of spectral line profiles in the H2CO and CS lines in our sample: double-peaked, 
with a peak and a "shoulder," and single symmetric lines. In addition, many sources have high 
velocity wings associated with unbound motions (outflows). 

The CS and II2CO line shapes differ from source to source, but are usually similar to each 
other. The N2H+ line, on the other hand, is Gaussian toward almost all sources. In sources with 
symmetric H2CO and CS lines, their peak velocity lies very close to that of the N2H^ line. In 
sources with double peaked H2CO and CS lines, the N2H+ peak velocity lies between the two 
peaks (or between the peak and the shoulder), indicating that the complex H2CO and CS line 
profiles arise from self-absorption at low velocities. As a proof of the lower optical depth of N2H"'" 
we derive optical depths r < 1 in the 101 — 012 line in 93% of the sources, and r < 1.6 in all of 
them. 

In most Class I sources, the H2CO and CS lines are broader than the N2H"'" line, but all three 
line profiles are essentially symmetric showing no evidence of self-absorption (Fig. ^). On the other 
hand, many of the observed CS and II2CO profiles toward the Class sources show evidence of 
self-absorption. The H2CO and/or CS spectral profiles are asymmetric with respect to the N2H+ 
center velocity in 15 out of 23 Class sources (Fig. [l|). Most of these asymmetric profiles have 
peaks whose velocity is significantly blue-shifted, rather than red-shifted, from the velocity of the 
optically thin N2H+ line. This high incidence of self-absorption toward Class sources is similar 
to what Gregersen et al. (1997) found in HCO^ lines. 

Besides the variation in line shapes discussed above, the line widths vary by a factor of five 
over our source sample. For example, the FWHM of the N2H"'" line varies from 0.3 km s^^ in 
Taurus to 1.8 km s^^ in Serpens. This is important to bear in mind when we use the peak 
velocity shift between lines to compare different sources in the following section. 

3.2. Quantifying Spectral Line Asymmetries 

To study the distribution of line profiles statistically, we need to quantify the observed line 
asymmetries. To do this, we define the non-dimensional parameter 



6V = {Vthick — Vthin) / ^Vthin, 



(1) 
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the velocity difference between tlie peaks of the optically thick and thin lines, normalized by the 
FWHM of the thin line. We use the normalized velocity difference 6V rather than the velocity 
difference Vthin ~ ^thick to reduce bias arising from lines having different width from one source 
to the next. This is a significant effect in our sample, as discussed above, and the normalization 
ensures that two lines of the same shape will have the same measure of asymmetry even if they 
have significantly different widths. To compute 5V , we first fit a Gaussian to the brightest peak 
of the optically thick line to define its velocity iVthick)- We then fit seven Gaussian hyperfine 
components to the N2H"'" J = 1 — line using the line parameters derived by Caselli et al. (1995); 
this gives the N2H+ velocity and FWHM {Vthin and /^thin)- 

We believe that 5V is a more convenient and robust measure of the line asymmetry than the 
even-odd component decomposition (Mardones et al. 1994), the shift in line centroid (Adelson, &; 
Leung 1988; Walker, Narayanan, & Boss 1994), the ratio of blue to red peak brightness (Myers 
et al. 1995), and the line skewness (Gregersen et al. 1997). The even-odd decomposition is too 
sensitive to the choice of reference velocity, and the line centroid is too sensitive to the choice of 
velocity window. The ratio of blue to red peak brightness cannot be used for lines where a second 
peak is not well defined. The skewness is much more sensitive than is 5V to asymmetry in the line 
wings, which may be due to outflow, not infall motions. On the other hand, the adopted 5V has 
the disadvantage that it is ambiguous for double-peaked lines with nearly equal peak intensities. 
In our sample this occurs in only 3 sources in the CS line (IRAS 16293-2422, L483 and Serp-S68N) 
and one source in the H2CO line (L483; see fig. ||). In the CS spectra, the ratio of the blue-red 
temperature difference to the rms noise ((T5 — Tr)/a) is 3.7, 3.0, and 0.8 in the sources IRAS 
16293-2422, L483, and Serp-S68N respectively, and in the L483 H2CO spectrum the same ratio is 
2.4. We regard a ratio of 2.0 or less as insignificant, so we drop the CS line in Serp-S68N from the 
sample. 

We investigated alternatives to fitting Gaussians to obtain the velocity of peak emission for 
the optically thick II2CO and CS lines. We explored measuring the velocity of the peaks directly 
from the spectra, calculating the moments of the spectra, and finding the velocity of the peak of 
the cross-correlation of the optically thick and thin line profiles. For simple lines all methods agree 
to within the rms uncertainty in center velocity due to noise. We adopt the direct fit method 
because it appears slightly more robust than the cross-correlation or first moment methods for 
sources with significant wing emission. In the optically thin N2II"'" J=l-0 line, the results of fitting 
seven Gaussian hyperfine components and of fitting a single Gaussian to the isolated 101 — 012 
line are consistent to within 0.01 km s~^, but the error in the former is smaller. 

Table | lists our main survey data: the observed velocities in the ortho-H2CO 2i2 — In, 
CS 2 — 1, and N2II"'" 1 — lines, and the values of 5Vh2CO = {yH2CO — '^N2H+) / ^'^N2H+ ai^d 
^Vcs = i^cs — ^N2H+) / ^^N2H+ computed from these velocities, as in equation |^. Table |^ lists 
the observed velocities in the smaller sample of para-Il2C0 2o2 — loi and C^^S 2 — 1 lines, which 
we observed as consistency checks. The velocity error estimates from the fitting routines are 
usually smaller than the channel width (except in low signal to noise spectra). We quote either one 
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half the channel width or the Gaussian estimate as the error in our peak velocity determinations, 
whichever is bigger. The error in Vthin ~ ythick is usually dominated by the error in the optically 
thick line profile velocity. This error is typically ~0.03 km s~^ and is always smaller than 0.11 
km s~^. 



3.3. Distributions of Normalized Velocity Differences 

Figure |^ shows histograms of the distribution of SVH2C0 and SVcs- The IRAM observations 
were used in figure ^ whenever sources were observed with both the SEST and IRAM telescopes 
because of the better spectral and spatial resolution (ie. each source is counted once). The 
histograms are clearly skewed towards negative velocities. This can also be seen from the 
means and standard errors of the mean of the distributions {6VH2C0) = —0.15 ± 0.08, and 
(SVcs) = —0.16 lb 0.07, given in Table |. We now explore the statistical significance of this result. 

We use the Student's t-test for a sample with unknown spread to compare the 5V distributions 
with a zero-mean normal distribution (our null hypothesis), and derive the probability (p) of the 
null hypothesis. The probability of drawing the distribution of 5V from a zero mean normal 
parent distribution is only 5% (3%) for the H2CO (CS) sample. 

Another way to quantify the asymmetries in the distribution of 5V is to simply count the 
excess (E) number of sources with negative compared to positive 6V . We define the "blue 
excess" E as (A^_ — N^)/N, where N- is the number of sources with a "significant" blue-shift 
{5V < —0.25), is the number of sources with a "significant" red-shift {6V > 0.25), Nq is the 
number of sources with no significant shift (—0.25 < 6V < 0.25), and is the total number 
of sources. Defined in this way, a large excess (E) corresponds to a small probability (p). We 
chose the threshold value 0.25 in the definition of E at about five times the typical rms error in 
6V (Table |3|) to screen out random contributions. The excess of sources with negative 6V in the 
histograms in figure |^ is 21% and 29% for the H2CO and CS samples respectively (25% overall). 

To understand the origin of the asymmetry of the 5V distribution, we evaluate statistics for 
subsamples defined by whether Tf^oi corresponds to a Class or a Class I source (from here on, 
we refer to sources with T^oi < 70 K as Class 0, and to sources with 70 < T^oi < 200 K as Class I). 
Table ^ lists (for the full sample and these subsamples) the quantities N; N^,Nq,N^; the mean; p, 
and E for the distributions of 6VH2CO and SVcs- 

Table || shows that the sample sources with Tboi < 70 K is entirely responsible for the 
histogram skewness in figure ^. Thus, Class sources have a probability of being drawn from a 
zero-mean normal distribution of 4%(0.4%) for the H2CO (CS) sample, while Class I sources have 
a probability of 81% (99%) for the H2CO (CS) sample. The overall excess of 25% of sources with 
negative 6V can be decomposed in a 47% excess (39% for H2CO and 55% for CS) among the Class 
sources, and 2% excess (4% and 0% respectively) among the Class I sources. This is illustrated 
in Fig. which shows the 6V distributions as a function of T^^i for the H2CO and CS lines. The 
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fraction of blue-shifted sources (negative 6V) in figure ^ is clearly greater for class than for Class 

1 sources. 

To summarize, sources with T^oi < 200 K (Class and I) tend to have asymmetric profiles in 
both of the optically thick lines that we used to trace high densities (> 10^ cm~^). In this sample 
there are ~25% more sources with a significant blue-shift than with a significant red-shift. This 
statistically significant tendency arises entirely from sources with Tboi < 70 K. In this Class 
subsample the blue-shifted excess is ~50%, while in the Class I subsample there is no significant 
blue- or red-shifted excess. 

3.4. Uncertainties Due to Choice of Line and to Angular Resolution 

Our statistical results are significant as long as they are independent of the instruments 
used, or the lines observed. We compare among different lines to better understand the sources 
of uncertainty in the statistical results discussed above, in particular possible effects of chemistry 
and resolution. We compare the optically thick line profiles of ortho and para H2CO, and CS; and 
the optically thin line profiles of N2II"'" and C^^S. 

3.4-1- Molecular Tracers as a Source of Uncertainty 

The H2CO and CS observations give statistically similar results as shown in section we 
now compare the line profiles individually in each source. The measured peak velocities in the CS 

2 — 1 and H2CO 2i2 — In lines agree quite well as can be seen by inspection of Table |3[ Figure |5| 
plots 5VH2CO against 6Vcs - Of the 39 sources observed in both lines, 4 (L1448mm, Serp-FIRSl, 
L1262 and L1527) deviate significantly from the line of perfect correlation. For the remaining 35, 
the best fit line has a slope of 1.05 =b 0.09, an intercept of —0.09 ± 0.04 and a correlation coefficient 
of 0.90. The small differences in the statistics derived above (see Table |^ arise mostly from the 
different CS and H2CO lines toward the four discrepant sources mentioned above. A paired t-test 
yields a probability of 0.91 that the means of the CS and H2CO 6V distributions are equal, so we 
treat both distributions as being drawn from the same parent distribution. We note that the CS 
and H2CO velocity differences were normalized by the same N2H^ line width, so the correlation 
in figure H is independent of the normalization. 

We observed the C^'^S 2 — 1 line toward 11 sources with the SEST and 2 sources with the 
IRAM telescope (table We evaluate the effect of using N2H+ instead of C^^S to compute 6V 
with CS 2 — 1 as the thick line. The difference 6V(j34 5 — 5V]\f^jj+ has mean it standard error of the 
mean (s.e.m.) of 0.06 ± 0.05, smaller than the width of one bin in figure ^. Most of the variation 
in 6V comes from the different widths of C^^S and N2H^ in some sources. The C^^S lines are 
broader than the N2H^ lines by a factor of 1.3 it 0.1. The velocity difference has mean it s.e.m 
0.02 ± 0.03 km s"^ {Vc'^'^s - yN2H+)- Only Serp-FIRSl and L1251B have C^^S velocities differing 
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by more than 0.1 km s ^ from the N2H"'" velocity. 

Ortho- and para-H2C0 are expected to trace the same physical conditions, and in LTE they 
are expected to have an ortho/para abundance ratio equal to three. The observed ortho and 
para line profiles are very similar in most sources, indicating high optical depth in both lines. In 
particular, the velocity of the peak is essentially the same in both transitions. The change in SV 
for all 12 sources has mean it s.e.m. 0.09 it 0.05 (para - ortho), smaller than the width of one bin 
in figure ^. The ortho- and para-H2C0 lines have different shapes only toward the sources L483 
and VLA1623, where the ortho-H2CO line peak velocity is more blue-shifted with respect to the 
N2H"'" velocity than is the para-H2C0 line. 

3.4-2. Angular Resolution as a Source of Uncertainty 

We explored the effects of different beam sizes by comparing the IRAM H2CO line profiles 
with average profiles over an effective 37" beam (3x3 points with 10" spacing and 16" beam) in a 
sample of 13 sources where we have made maps (Mardones et al. 1997). The self-absorption depth 
and blue to red peak ratio vary somewhat, but the peak velocities are consistent to within 0.1 
km s~^. The only exception is Serpens-FIRSl, where the red peak is brighter in the IRAM central 
spectrum but the blue peak is brighter in the IRAM 3x3 average. The difference in 5V between 
the average spectrum of the 3x3 map and the center spectrum has mean it s.e.m. 0.04 it 0.03. In 
general, then, the measured 6V is not very sensitive to variations in pointing or spatial resolution. 

We also compared SEST and IRAM observations to study the variation of the line profiles 
when changing spectral and spatial resolution. The observed SEST line profiles are on average 40% 
fainter than the corresponding IRAM line profile. This is likely due to beam dilution indicating 
that the emitting gas is spatially concentrated on scales smaller than the SEST main beam (36"). 
Some of the sources proved to be very compact as judged by the high IRAM/SEST intensity ratio 
(notably 16293-2422 and B335). The shapes of the lines observed with both telescopes are very 
similar in most of the sources, with a few exceptions noted below. For sources with double-peaked 
H2CO line profiles the overall intensity and width of the lines is similar at both telescopes, but the 
self-absorption depth and the ratio of the two peaks can vary significantly. The dips are always 
deeper in the IRAM spectra, maybe indicating that the absorption is also diluted in the SEST 
beam. The peak velocities measured in the SEST and IRAM spectra are consistent with each 
other. The change in 5V using the II2CO velocities measured at the SEST and IRAM has mean 
zb s.e.m. 0.01 lb 0.05 (SEST-IRAM), for the 10 sources in common, excluding Serp-FIRSl. The 
SEST spectrum toward Serp-FIRSl is almost identical to the IRAM 3' x 3' average spectrum. We 
find no telescope biases in our observations that could affect significantly the statistics. 



- 12 - 



4. Discussion 

We have found an excess of sources with asymmetric optically thick lines (double-peaked or 
shoulder profiles) toward a sample of 47 nearby low mass YSOs. This asymmetry is such that the 
peak velocity of the optically thick line is preferentially blue-shifted with respect to the optically 
thin velocity. This blue-shifted excess is statistically more significant in subsamples of sources 
with Thoi < 70 K. The observed blue excess of dV values is similar for the CS and H2CO tracers 
and is independent of the telescope used. It is, therefore, an inherent property of the sources in 
the samples. 



4.1. Kinematic Interpretation of the Observed Velocity Differences 

The observed N2H+ line profile is single in all sources observed, ruling out multiple 
components as the cause of the asymmetric lines observed. In a few cases (eg. L1448mm), the 
N2H+ line profile departs significantly from a Gaussian. However, in those cases the peak of the 
N2H"'" line profile lies within the H2CO or CS absorption dip and does not coincide with the peak 
of the H2CO or CS profile. 

A foreground absorbing layer could reproduce the asymmetric line profiles we observe. 
However, we do not expect such a layer of dense gas (having the systemic velocity of the source) 
to lie in front of most of the sources observed unless it is physically associated. On the other hand, 
the spatial variation of optically thick line profiles show that the asymmetries are concentrated 
toward the sources (§ |3.4.2 ; also Zhou et al. 1993, Myers et al. 1995, Mardones et al. 1997) making 



the casual superposition of two clouds even more unlikely. We are led to explain the variety of line 
profiles observed as arising from local self-absorption, i.e. absorption by physically associated gas. 

We now ask what kinematics should the local absorbing gas have with respect to the YSO to 
reproduce the observed 6V distributions. If the kinematics around a YSO are dominated by infall, 
then we expect the foreground gas to be preferentially red-shifted with respect to the systemic 
velocity of the source. If this foreground gas has colder excitation temperature than the dense 
gas in the inner cloud, it will absorb preferentially at red-shifted velocities causing asymmetric 
line profiles (eg Leung &i Brown 1977; Zhou 1992). Thus it will have a distribution of 5V skewed 
toward blue velocities as observed. If systematic motions are due solely to infall toward a high 
fraction of Class sources, these motions can explain the observed statistics. 

The observed blue excess of 6V can arise from infall, as opposed to expansion; but other 
motions, including rotation and bipolar outflows, need to be considered. Adelson & Leung (1988) 
showed that rotation can cause infall-like and expansion-like line profiles on either side of the 
rotation axis (see also Zhou 1995). The observed effect depends on the orientation of the rotation 
axis with respect to the line of sight. A source sample with randomly oriented rotation axes 
would produce a symmetric distribution of 6V. Similarly, bipolar outflows also cause infall- and 
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expansion-like line profiles on either side of the source (Cabrit, & Bertout 1986). An ensemble 
of bipolar outflows with random orientation would also produce a symmetric distribution of 6V, 
unlike our results. Thus our statistical results do not rule out the contribution of outflow and 
rotation to the line profiles, but they require inward motions in addition to any outflow and 
rotation in a significant fraction of our source sample. 

Outflow motions greatly affect the observed line profiles toward some sources, in particular 
among the Class sample. If outflow axes in this sample were not distributed randomly, or if the 
sample were too small, the outflows could, in principle, produce the observed statistical effect. We 
now examine on a source-by-source basis whether the observed outflow wings are associated with 
the excess of sources with blue 5V in our sample. We subtracted one or two best-fit Gaussians 
from the H2CO and CS line peaks to measure the remaining blue and red wings in every spectrum. 
We then labeled each source's wings as blue if the blue wing is >25% brighter than the red, red 
if the red wing is >25% brighter than the blue, and nil otherwise. We find as many sources with 
brighter blue and red wings (~ 25% of each) in both the H2CO and CS lines, therefore the excess 
of sources with negative 6V is not associated with a similar excess of sources with bright blue 
wings. 

To test this further, we defined a wing parameter with values -1, 0, or 1 for sources with 
brighter wings to the blue, neither, or red respectively, and a line-core parameter with values -1, 0, 
or 1 for sources with a measured 6V < —0.25, —0.25 < 5V < 0.25, or 5V > 0.25 respectively. We 
find a Spearman rank correlation coefficient between the wing and the line-core parameters of 0.0 
for the H2CO sample and -0.1 for the CS sample. We conclude that the observed wing emission is 
not correlated with, and thus is not responsible for the excess of sources with negative 6V derived 
in section 3. 

Finally, it is possible that star forming groups could have large scale rotation or magnetic 
fields imposing a characteristic orientation to the disks and outflows around YSOs. This effect 
could bias the statistics within each cluster. This is evidently not the case in this survey, because 
the projected outflow axes have a wide distribution in position angle within each cluster. The 
distribution of axis position angles of all outflows with well known orientation in NGC1333 has 
an rms of 56° (Bally, Devine, & Reipurth 1996). Bontemps et al. (1996) found a similar range of 
outflow axis position angles, with an rms of 55°in Ophiuchus. In Serpens, the position angles of 
the outflows axis of the sources S68N, FIRSl, SMM2, SMM4, and SMM5 in Serpens have an rms 
of 43°(White, Casali,, & Eiroa 1995, Williams et al. 1997). Therefore, within each of the clusters 
in our survey we can regard the source outflow axes as randomly oriented. 



4.2. Infall Toward Class and Class I Sources 



We have shown that the observed distribution of 6V is explained best by motions which 
include infall in addition to outflow and/or rotation. Furthermore, we showed in §3.5 that in our 
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sample the infall signatures are seen almost entirely toward Class rather than Class I sources. 
The most straightforward interpretation of these results is that Class sources differ from Class 
I sources in having a much greater incidence of inward motions for the conditions which our 
observations probe -size scales of a few x 0.01 pc and gas density of about 10^ cm~^. This 
difference, which has not been demonstrated before, adds support to the idea that Class sources 
have not yet accreted most of their mass, based mainly on the greater submillimeter continuum 
emission of Class than Class I sources (Andre et al. 1993). 

The lack of infall signatures in Class I objects, however, does not necessarily mean that 
infall is absent in these sources. Adams, Lada, k, Shu (1987) showed that to model the spectral 
energy distribution of Class I sources it is necessary to include the contribution of an infalling 
envelope. The lack of blue-shifted line profiles toward Class I sources could therefore result from a 
lower sensitivity of our observations to the physical conditions in the infalling envelopes of Class I 
sources. This would occur if the circumstellar envelopes in Class I sources have less optical depth 
than in Class sources, or if the circumstellar envelopes in Class I sources have smaller excitation 
temperature gradients than in Class sources (in the optically thick CS and H2CO lines). 

To test the above ideas, we examined the optical depths determined from the multicomponent 
fits to the N2H+ spectra for the Class and Class I sources in our sample. We found no significant 
difference between the two groups: the N2H"'" optical depth, summed over all seven hyperfine 
components, has mean and s.e.m. 4.9 it 0.6 for the Class 0, and 4.6 it 0.7 for the Class I source 
sample. This suggests that the greater incidence of blue-shifted profiles in Class than in Class I 
sources does not arise from differences in line optical depth between the two samples. However, 
the lack of self-absorption features in the Class I source line profiles may indicate that the CS and 
H2CO optical depths could still be higher in Class than in Class I sources. Also, there are more 
weak single-peaked sources in the Class I than in the Class source samples. 

Sources whose outflow axes are in the plane of the sky, having edge on disks, may have 
signiflcant obscuration to optical and near-infrared radiation, lowering their observed Thoi- This 
effect would tend to mix the Class and Class I sources. However, the circumstellar envelopes 
are most likely optically thin at sub-mm wavelengths, so a sub-mm wavelength excess should 
distinguish Class from Class I sources, independent of orientation. Among our sample of Class 
sources, only 03282-1-3035 and L673A have not been detected at sub-mm wavelengths; moving 
them to the Class I sample would only emphasize the observed difference between the Class and 
Class I samples. 

Thus, we found a significant difference in the incidence of blue-shifted line asymmetry between 
the Class and Class I sources, indicating different physical conditions in the circumstellar 
envelopes around Class and Class I sources. However, our observations do not necessarily rule 
out the presence of inward motions onto Class I sources; tracers with higher optical depths may 
be needed to probe the circumstellar kinematics around Class I sources. 
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4.3. Identification of Infall Candidates 

We summarize in table ^ the evidence for infall, expansion, or neither for the sources in our 
sample. We list all 23 sources with \SV\ > 0.25 in either the H2CO or CS lines. We label the lines 
b if they have blue asymmetry {6V < —0.25), r if they have red asymmetry {6V > 0.25), or n if 
they have neither (-0.25 < 5V < 0.25). We include IRAS 13036-7644 in the CS sample, using 
the spectrum published by Lehtinen (1997). The correspondance between both lines is extremely 
good, as expected from figure |5[ Following Gregersen et al. (1997), we identify infall candidates 
as those sources with evidence for infall (at least one "b") and no evidence for expansion in any 
line (no "r"s). Therefore, from the H2CO and CS observations presented in this paper, we select 
15 sources as infaU candidates: IRAS 03256+3055, NGC1333-4A, NGC1333-4B, L1527, IRAS 
13036-7644, VLA1623, WL22, IRAS 16293-2422, L483, S68N, SMM5, SMM4, B335, L1157, and 
L1251B. These are Class sources {Thoi < 70 K) with the exception of IRAS 03256+3055, WL22 
and L1251B. In contrast to the 15 (CS and II2CO) sources with at least one b and no rs, there are 
only 4 with the opposite property -at least one r and no bs (NGC1333-2, L1551-5, L1551NE and 
L43). 

Among our infall candidates, L1527, IRAS 16293-2422 and B335 have been widely studied 
(eg. Walker et al. 1986; Zhou et al. 1993, 1994; Mardones et al. 1994; Myers et al. 1995; Ohashi 
et al. 1997; Gregersen et al. 1997). They represent the best cases for infall to date due to the 
agreement between kinematic models and observations. The sources NGC1333-2 (Ward- Thompson 
et al. 1996), IRAS 13036-7644 (Lehtinen 1997), L483 (Myers et al. 1995), S68N (Hurt, Barsony, k 
Wootten 1996), SMM4 (Hurt et al. 1996, Gregersen et al. 1997), L1157 (Gueth et al. 1997), and 
L1251B (Myers et al. 1996) have been recently identified as infall candidates but have not yet been 
studied as thoroughly as L1527, IRAS 16293-2422, and B335. The sources IRAS 03256+3055, 
NGC1333-4A and 4B, VLA1623, WL22, and SMM5 are new kinematic infall candidates based on 
this work. 

4.4. Comparison with Previous Surveys 

Wang et al. (1995) observed 12 globules at the CSO using the H2CO 3i2 — 2ii hue looking 
for kinematic infall signatures. They do not find convincing spectral evidence of infall toward any 
source. This may be because none of their sources are known to be Class YSOs; or due to the 
expected lower optical depth of the H2CO 3i2 — 2ii line. 

Gregersen et al. (1997) observed 23 Class sources in the HCO^ 3 — 2 and 4 — 3 lines looking 
for infall signatures. Their sample contains sources more distant than 400 pc, but the overlap with 
our sample is still considerable, with 15 sources in common. We derived SV from the HCO"*" and 
H^^CO''" spectra presented by Gregersen et al. toward 13 of those sources, and assigned asymmetry 
labels b,r, and n to them in the same way as for the CS and H2CO spectra. In table |6| we list the 
asymmetry labels from available observations: from this paper, from Gregersen et al. (1997), from 
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Ward-Thompson et al. (1996), and from Lehtinen (1997). The correspondance among all lines is 
generally good. Seven of the nine sources with infall asymmetry in the HCO^ lines have the same 
sense of asymmetry in the H2CO and CS lines. On the other hand, the two sources with expansion 
asymmetry in the HCO"*" lines show infall asymmetry in the H2CO and/or CS lines (L1448mm 
and L483), and one source (NGC1333-2) shows infall asymmetry in HCO^ 4 — 3 (Ward-Thompson 
et al 1996) and expansion asymmetry in CS and H2CO. These disagreements may reflect differing 
degrees of outflow contamination from one line to the next. 

Assuming as in |4.3| that the best infall candidates are those sources with evidence for infall 
in at least one line, and with no evidence for expansion in any line, we consider again the sources 
in Table now according to their asymmetry in HCO^, as well as in CS, H2CO. Adding this 
HCO"*" data and applying this procedure removes only one source, L483, from the CS and H2CO 
list given in §^^, leaving the 14 sources indicated in Table ^ as the best infall candidates. These 
sources have substantial internal agreement among their asymmetry labels: 5 have 4 bs, 2 have 3 
bs, 5 have 2 bs, and 2 have 1 b. At present, these sources have the most evidence from available 
spectral line surveys for infall, and they have the least evidence against infall. 



The kinematic infall candidates identified in this paper include two sources -B335 and L1527- 
which exemplify in several ways the "standard model" of low-mass star formation (e.g. Shu, 
Adams, &: Lizano 1987). These sources are relatively isolated and have NH3 line widths dominated 
by thermal, rather than turbulent motions (Myers, & Benson 1983; Benson, & Myers 1989). Their 
profiles of H2CO and other lines are well- fit by radiative transfer models whose kinematics are 
specified by the "inside-out" mode of collapse of a singular isothermal sphere (Shu 1977). 

But these two sources are not representative of the properties of the 15 kinematic infall 
candidates identified by the CS and H2CO lines in this paper. Most of these sources differ from 
the standard model in two ways: their line widths are dominated by turbulent motions, and many 
of them are found in stellar groups rather than in isolation. 

We quantify the turbulent motions of the infall candidates by comparing the FWHM of the 
optically thin N2H+ line, given in Table ^, with that expected for gas having equal thermal and 
nonthermal motions. 



where T is the gas kinetic temperature, mobs is the mass of the observed species (29 amu for 
N2H"'"), and (m) is the mean molecular mass (2.3 amu). For gas with T = 15 K, 9 out of 15 
sources, or 60%, have N2H^ line widths AV > AVq = 0.57 km s^^, i.e. having greater non-thermal 
than thermal motions. If T = 10 K, this fraction increases to 12/15, or 80%. Clearly most of the 
kinematic infall candidates have non-thermal motions exceeding their thermal motions. 



4.5. 



Thermal and Turbulent Infall 




(2) 
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The non-thermal motions in the N2H"'" hne width could arise from motions which appear 
random to our resolution, which we designate as turbulence, and also from systematic motions, 
including outflow, infall, and rotation. Since N2H"'" lines (like NH3 lines) tend to trace quiescent 
core structure (Womack et al. 1992) and do not track outflow maps (Caselli et al. 1997), we expect 
that the outflow contribution to the N2H"'" line width is relatively small compared to the turbulent 
component. Similarly, the relative contribution of rotation and infall motions to the optically thin 
line width is significantly smaller than the random contribution, according to infall models which 
match the profiles of infall candidates (e.g. Zhou et al. 1993). 

To assess the proportion of the infall candidates associated with groups or clusters, we 
consider those infall sources in the NGC1333, Ophiuchus, and Serpens clouds, which are associated 
with a substantial cluster of embedded young stellar objects. Of the 15 infall candidates listed in 
Section § |4.3| , 8 or 53% belong to these three clusters. This accounting is conservative in that we 
do not count L1251B which belongs to a smaller group of at least 5 near infrared sources (Hodapp 
1994). Of these 8 sources in clusters, most (6) have N2H+ line widths dominated by turbulent 
motions. 

This relatively high incidence of turbulent motions and cluster sources in our list of kinematic 
infall candidates implies that the standard model of isolated star formation in a thermally 
dominated core is not representative of most of our infall candidates. This is an area where more 
theoretical models are needed. Models of star formation which account for significant turbulent 
motions as part of their initial conditions have been presented for a spherically symmetric 
geometry and for the formation of a single star (Caselli &; Myers 1995; McLaughlin & Pudritz 
1997); but these models lack a physical basis for the non-thermal motions, and do not account for 
the formation of multiple stars. 



4.6. Future Prospects 

The infall candidates identified in this paper are good candidates for the study of star-forming 
infall. However, it is not established from the survey data presented here that their motions are 
necessarily either "star- forming" or "gravitational." 

A further step to establish if the blue-shift identified in Class sources arise from infalling 
envelopes is to map them: the blue-shifted, asymmetrical line profiles should be spatially 
concentrated toward the source. Furthermore, maps are needed to disentangle the effects of the 
outflows. Maps of observed line asymmetry (eg. Adelson & Leung 1988; Walker et al. 1994) can 
distinguish whether infall dominates (at low velocities) over rotation or expansion in a given 
region. 

To establish whether the observed infall kinematics are gravitational, the observations should 
be compared with kinematic models, which predict the spatial variation of the line proflles (eg. 
Zhou et al 1993). Furthermore, continuum observations may help determine the YSO and envelope 
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masses to compare with model predictions. High spatial resolution attainable with interferometers 
can locate the infall regions in relation to the position of the continuum peak of the YSO. All 
of these together should constrain significantly kinematic models of gravitational collapse, to 
determine better the physical basis of the inward motions inferred in this paper. In addition, 
further models and observations may be needed to resolve the apparent conflicts presented in table 
^. Such models will have to incorporate both infall and outflow kinematics. 

5. Conclusions 

Among 47 of the youngest nearby low mass PMS stars (spectral classes 0, or T^oi < 70 K, and 
I, or 70 < Tboi < 200 K) a significant fraction show complex line profiles with double peaks or a 
peak and a shoulder in the optically thick CS 2 — 1 and H2CO 2i2 — In transitions, which trace 
dense gas (n > 10"^ cm~^). 

These complex spectral profiles are not due to multiple components along the line of sight, 
since the optically thin N2H+ 101 — 012 line, observed toward the same sources, is single. The 
N2H+ line velocities are consistent with the C^^S 2-1 velocities. 

In both CS or II2CO samples, about half the sources have a significant velocity-shift between 
the optically thick and thin lines. 

In both CS or II2CO samples, there are ~25% more sources whose velocity shift is to the 
blue than to the red. This statistically significant tendency arises entirely from sources with 
TfeoZ < 70 K (Class sources). In this Class sample the blue-shifted excess is 50% (9 out of 23 
for II2CO and 10 out of 19 for CS), while in the Class I subsample there is no significant blue- or 
red-shifted excess. 

The observed excess of blue-shifts cannot be reproduced by a sample of sources whose 
kinematics are dominated either by bipolar outflows or rotation, if their symmetry axes are 
oriented in random directions. On the other hand, infall motions alone, or infall combined with 
rotation and/or outflows, can explain the excess of blue-shifts. 

We identify a total of 15 spectroscopic infall candidates, of which 6 have not been identified 
previously. Most of these candidates differ from the well-known sources B335 and L1527, which 
exemplify the "standard model" of an isolated core whose line width is dominated by thermal 
motions. They tend to have primarily turbulent line-widths, or are associated with embedded 
young clusters, or both. 
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TABLE 1 
Source List 



Source^ 


Name 


RA (1950.0) 

K m. s 


DEC (1950.0) 

O ( (( 


Distance 
pc 


Ti,„, 
K 


Na 


References'^ 


Outflow' 


03225+3034 


L1448-IRS3 


3 22 31.5 


30 34 49 


350 


70 


5 


1,2 


CO 


L1448mm 




3 22 34.4 


30 33 35 


350 


56 


6 


3 


CO,HH 


03235+3004 




3 23 33.0 


30 04 59 


350 


136 


5 


1,2 




03245+3002 


L1455, RN015 


3 24 34.9 


30 02 36 


350 


71 


5 


1,2 


CO 


03256+3055 




3 25 39.2 


30 55 20 


350 


74 


4 


4 




03259+3105 


SS V13 


3 25 57.9 


31 05 50 


350 


90 


6 


2,5 


CO,HH 


NGC1333-2 




3 25 49.9 


31 04 16 


350 


52 


9 


6,7 


CO,HH 


NGC1333-4A 




3 26 04.8 


31 03 13 


350 


34 


8 


6,8 


CO 


NGC1333-4B 




3 26 06.5 


31 02 51 


350 


36 


8 


6,8 


CO 


03282+3035 




3 28 15.2 


30 35 14 


350 


35'' 


2 


1 


CO,HH 


HH211mm 




3 40 48.7 


31 51 24 


350 


30*^ 





9 


CO,HH 


04166 + 2706 


B213NW 


4 16 36.0 


27 06 00 


140 


139 


8 


10 


CO 


04169 + 2702 


B213NW 


4 16 54.0 


27 01 59 


140 


170 


11 


10 


CO 


04287+1801 


L1551-IRS5 


4 28 40.2 


18 01 42 


140 


95 


19 


10 


CO,HH 


04289 + 1802 


L1551NE 


4 28 50.8 


18 02 11 


140 


75 


8 


10 




04325 + 2402 


L1535 


4 32 33.5 


24 02 15 


140 


157 


5 


11 


CO 


04365 + 2535 


TMCIA 


4 36 31.2 


25 35 56 


140 


170 


11 


10 


CO 


04368+2557 


L1527 


4 36 49.3 


25 57 16 


140 


60 


11 


10 


CO,HH 


11051-7706 


CedllO-IRS4 


11 05 23.1 


-77 06 14 


150 


95 


9 


11 


CO,HH 


11590-6452 


Sal36, BHR71 


11 59 03.1 


-64 52 11 


150 


56 


5 


1,12 


CO,HH 


12553-7651 




12 55 20.3 


-76 51 22 


150 


99 


4 


11 




13036-7644 


Sal60, BHR86 


13 03 41.4 


-76 44 03 


150 


63 


5 


11 


CO 


15398-3359 


B228 


15 39 50.4 


-33 59 42 


150 


61 


4 


11 


HH 


15420-3408 


Sz68 


15 42 01.3 


-34 08 09 


150 


120 


5 


1,5 


HH 


VLA1623 




16 23 24.9 


-24 17 46 


160 


30'' 


4 


13 


CO,HH 


WL22 


IRS27, YWL4 


16 23 57.3 


-24 28 15 


160 


86 


10 


10 




16244-2432 


IRS44, YLW16 


16 24 26.2 


-24 32 53 


160 


157 


10 


10 


CO 


16293-2422 




16 29 20.9 


-24 22 16 


160 


43 


6 


10 


CO 


16316-1540 


L43, RN091 


16 31 37.7 


-15 40 52 


160 


79 


6 


1,14 


CO 


16544-1604 


L146, CB68 


16 54 27.2 


-16 04 48 


160 


74 


3 


1 


CO 


18148-0440 


L483 


18 14 50.6 


-04 40 49 


200 


50 


8 


1,14,15 


CO,HH 


Serp S68N 




18 27 15.2 


01 14 57 


310 


40'' 


4 


16,17 


CO 


Serp FIRSl 




18 27 17.4 


01 13 16 


310 


51 


7 


16,17 


CO 


Serp SMM5 




18 27 18.9 


01 14 36 


310 


59 


7 


16,17 


HH 


Serp SMM4 




18 27 24.3 


01 11 11 


310 


43 


6 


16,17 


CO 


Serp SMM3 




18 27 27.3 


01 11 55 


310 


46 


5 


16,17 


CO 


18331-0035 


HH108/109IRS 


18 33 07.6 


-00 35 48 


310 


54 


4 


1,5 


HH 


CrA-IRAS32 




18 59 35.8 


-37 11 53 


130 


148 


4 


11 




19156 + 1906 


L723 


19 15 41.3 


19 06 47 


200 


50 


4 


1,5 


CO,HH 


19180 + 1114 


L673A 


19 18 04.6 


11 14 12 


300 


55 


4 


1,14 


CO 


19345 + 0727 


B335, CB199 


19 34 35.7 


07 27 20 


250 


37 


10 


1,18 


CO,HH 


20353+6742 


L1152 


20 35 19.4 


67 42 30 


350 


72 


3 


1 


CO 


20386 + 6751 


L1157 


20 38 39.6 


67 51 33 


350 


62 


3 


1 


CO,HH 


21017+6742 


L1172 


21 01 44.2 


67 42 24 


350 


72 


5 


1,14 


CO 


22343+7501 


L1251A 


22 34 22.0 


75 01 32 


300 


108 


7 


1 


CO,HH 


22376 + 7455 


L1251B 


22 37 40.8 


74 55 50 


300 


91 


4 


1 


CO,HH 


23238+7401 


L1262, CB244 


23 23 48.7 


74 01 08 


250 


104 


7 


1,19 


CO 



^IRAS name whenever applicable. 

'^References to SED observations (or T^oO only. 

'CO indicates that an outflow has been detected as high velocity molecular gas (White et al. 1994; Wu et al. 1996; Bontemps et al. 1996; Williams 
et al. 1997). HH indicates that there is evidence of shocked molecular gas from optical or near-infrared line emission (Bally et al. 1993a, 1993b; Bourke 
1994; Reipurth 1994; Dent et al. 1995; Fuller et al. 1995; Devine et al. 1997; Eiroa et al. 1997) 

''educated guess, we estimate limits Ti,„, < 50 K (03282+3035), Ti,„! > 8 K (HH211 & VLA1623) and Ti,„, < 80 K (Serp 68N) 

REFERENCES. — (1) IRAS Point Source Catalogue; (2) Ladd et al 1993; (3) Bachiller et al 1991; (4) Clark 1991; (5) Reipurth et al 1993; (6) Jennings 
et al 1987; (7) Sandell et al 1994; (8) Sandell et al 1991; (9) McCaughrean et al 1994; (10) Chen et al 1995; (11) Chen et al. 1997; (12) Bourke et al 
1995; (13) Andre et al 1993; (14) Ladd et al 1991; (15) Fuller et al 1995; (16) Casali et al 1993; (17) Hurt & Barsony 1996; (18) Chandler et al 1990; 
(19) Yun & Clemens 1995 
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TABLE 2 

Sources not detected in either the H2CO line, the N2H+ line, or in both lines 



Source^ 


Name 


RA (1950.0) 

h jn s 


DEC (1950.0) 


Distance 


Ti,t>i 


References^ 


Line(s) 


rms noise^ 


Telescope 






/ // 


pc 


K 






K 




03382+3145 




3 38 15.0 


31 45 11 


350 


85 


1,2 


H2CO 




0.15 


IRAM 


03385+3149 




3 38 30.9 


31 49 24 


350 


110 


1,2 


H2CO, N2H+ 





.15, 0.10 


IRAM 


04108+2803 


L1495 


4 10 49.3 


28 03 00 


140 


205 


3 


H2CO, N2H+ 





.15, 0.11 


IRAM 


04302 + 2247 


L1536 


4 30 16.4 


22 47 04 


140 


202 


3 


H2CO, N2H+ 





.13, 0.10 


IRAM 


04361 + 2547 




4 36 06.6 


25 47 00 


140 


144 


3 


N2H+ 




0.11 


IRAM 


Cham B35 




11 05 55.0 


-77 06 18 


150 


57 


4 


N2H+ 




0.04 


SEST 


Cham C9-2 




11 07 12.2 


-77 27 37 


150 


175 


4 


N2H+ 




0.05 


SEST 


Cham IRN 




11 07 14.1 


-77 27 28 


150 


208 


4 


H2CO, N2H+ 





.12, 0.05 


SEST 


11306-6311 


DC294. 3-2.0 


11 30 42.5 


-63 10 45 


150 


71 


4 


H2CO, N2H+ 





.12, 0.04 


SEST 


11436-6017 




11 43 35.2 


-60 17 27 


150 


103 


4 


N2H+ 




0.05 


SEST 


12500-7658 




12 50 03.3 


-76 58 39 


150 


94 


4 


H2CO, N2H+ 





.10, 0.04 


SEST 


12533-7632 




12 53 19.3 


-76 32 12 


150 


95 


4 


H2CO, N2H+ 





,09, 0.04 


SEST 


15215-6056 


DC320.5-3.6 


15 21 30.7 


-60 56 19 


150 


84 


1,5 


N2H+ 




0.03 


SEST 


16017-3936 




16 01 42.5 


-39 36 56 


150 


99 


4 


H2CO, N2H+ 





.05, 0.04 


SEST 


GSS30-IRS1 




16 23 19.9 


-24 16 18 


160 


175 


3 


N2H+ 




0.10 


SEST 


Oph IRS8 




16 22 46.3 


-24 11 43 


160 


203 


3 


H2CO, N2H+ 





.15, 0.13 


IRAM 


16245-3859 




16 24 30.2 


-38 59 07 


150 


129 


4 


H2CO, N2H+ 





.04, 0.04 


SEST 


17151-3202 


DC354.2+3.2 


17 15 08.6 


-32 02 24 


170 


134 


1,5 


H2CO, N2H+ 





.05, 0.04 


SEST 


19180 + 1116 


L673B 


19 18 01.3 


11 16 27 


300 


94 


1,6 


H2CO, N2H+ 





.17, 0.19 


IRAM 



^IRAS name whenever applicable 

^References to SED observations (or Tt,o;) only 

^noise in units of T^t- 

REFERENCES. — (1) IRAS Point Source Catalogue; (2) Ladd et al 1993; (3) Chen et al 1995; (4) Chen et ah 1997; (5) Persi et al 1990; (6) Ladd et 
al 1991 
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TABLE 3 





H2CO 2 


12 — 111, 


CS 2 


— 1 AND 


N2H + 


1 - 


Spectral 


Line Velocities 










source 


Telescope 


y 1^2 CO 




Vcs 








^yN2H+ 














km s 


1 




km s~ ^ 




km s~ ^ 






km s ^ 














L1448-IRS3 


IRAM 


4.68 ± 





.02'' 


4,70 + 


,02 




52 + 


01 





96 + 


1 





,16 + 





,02 





,19 + 0,02 


L1448mm 


IRAM 


5.59 ± 





.03 


4,39 + 


,02 




02 + 


02 





84 + 


05 





,67 + 





,06 


-0 


,75 + 0,06 


03235+3004 


IRAM 


5.18 + 





.02 


5,19 + 


,03 


5 


,14 + 


01 





,35 + 


, 1 





,11 + 





,06 





,15 + 0,09 


L1455 


IRAM 


4.72 ± 





.02 








78 + 


01 





69 + 


1 


-0 


,09 + 





,03 






03256+3055 


IRAM 


7.32 ± 





.02 


7,28 + 


,02 




39 + 


01 





34 + 


1 


-0 


,21 + 





,07 


-0 


,31 + 0,07 


SSV13 


IRAM 


8.55 ± 





.02 


8,36 + 


,03 







01 





v> -t- n 


n 1 





,14 + 





,03 


-0 


,13 + 0,04 


NGC1333-2 


IRAM 


7.98 ± 





.03 


7,71 + 


,02 


1 


,52 + 


01 





, 94 + 


, 1 





,49 + 





,03 





,20 + 0,02 


NGC1333-4A 


IRAM 


6.63 ± 





.02 


6,73 + 


,02 


1 


,25 + 


01 


1 


,32 + 


,02 


-0 


,47 + 





,02 


-0 


,39 + 0,02 


NGC1333-4B 


IRAM 


6.60 ± 





.02 


6,75 + 


,03 




1 1 -t- n 


01 


1 




n'> 

, 


-0 


,49 + 





,02 


-0 


,34 + 0,03 


03282+3035 


IRAM 


7.02 ± 





.02 


6,99 + 


,02 


Q 


96 + 


01 





47 + 


03 





,13 + 





,05 





,06 + 0,05 


HH211MM 


IRAM 


9.16 ± 





.02 


9,07 + 


,02 


g 


10 + 


01 





45 + 


02 





,13 + 





,05 


-0 


,06 + 0,05 


04166 + 2706 


IRAM 


6.62 ± 





.02 


6,65 + 


,02 


D 


, Do ^ u 


01 







n 1 


-0 


,16 + 





,07 


-0 


,08 + 0,07 


04169 + 2702 


IRAM 


6.90 ± 





.02 


6,88 + 


,02 


Q 


91 + 


01 





29 + 


1 





,03 + 





,08 


-0 


,11 + 0,08 


L1551-IRS5 


IRAM 


6.39 ± 





.02 


6,77 + 


,02 


Q 


42 + 


01 





80 + 


1 


-0 


,04 + 





,03 





,44 + 0,03 


L1551NE 


IRAM 


7.02 ± 





.02 


7,05 + 


,02 


Q 


65 + 


,01 





47 + 


03 





,78 + 





,07 





,85 + 0,07 


04325 + 2402 


IRAM 


5.65 ± 





.02 


5,73 + 


,02 




67 + 


,01 





40 + 


1 


-0 


,06 + 





,06 





,14 + 0,06 


04365 + 2535 


IRAM 


6.45 ± 





.04 


6,46 + 


,02 


Q 


40 + 


01 





42 + 


1 





,11 + 





,05 





,14 + 0,05 


L1527 


IRAM 


5.59 ± 





.02 


5,91 + 


,02 


c 



Q 1 -t- n 


,01 







n 1 


-0 


,99 + 





,08 





,01 + 0,07 


11051-7706 


SEST 


4.20 ± 





.11 








33 + 


01 





63 + 


1 


-0 


,20 + 





,18 






11590-6452 


SEST 


-4.59 ± 





.02 






^ 


46 + 


,01 





88 + 


03 


-0 


,15 + 





,03 






12553-7651 


SEST 


3.24 + 





.05 








22 + 


,01 





46 + 


1 





,05 + 





,11 






13036-7644 


SEST 


3.36 ± 





.03 






3 


, 76 + 


01 


1 


,01+0 


,02 


-0 


,40 + 





,03 






15398-3359 


SEST 


5.05 ± 





.02 








14 + 


01 





49 + 


1 


-0 


,17 + 





,05 






15420-3408 


SEST 


5.09 ± 





.02 








04 + 


,01 





58 + 


03 





,10 + 





,04 






VLA1623 


IRAM 


3.31 ± 





.02 


3,27 + 


,02 




67 + 


01 





56 + 


1 


-0 


,64 + 





,04 


-0 


,71 + 0,04 




SEST 


3.15 ± 





.02 








65 + 


,01 





71 + 


1 














WL22 


IRAM 


3.50 ± 





.02 


3,66 + 


,02 


Q 


,00 ^ u 


,01 





,00 ^ u 


n 1 


-0 


,99 + 





,07 


-0 


,53 + 0,07 




SEST 


3.59 ± 





.02 






3 


80 + 


01 





53 + 


1 














16244-2432 


SEST 


4.28 ± 





.02 








14 + 


01 





72 + 


05 





,19 + 





,03 








IRAM 


3.12 ± 





.02 


3,25 + 


,02 




16 + 


,01 





78 + 


1 


-1 


,34 + 





,03 


-1 


,16 + 0,03 


16293-2422 


SEST 


3.09 ± 





,02 






4 


,04+0 


01 


1 


, 13 + 


,02 














L43 


IRAM 


0.46 ± 





.02 


0,64 + 


,02 





,44 + 


01 





,48 + 


,01 





,04 + 





,05 





,40 + 0,05 




SEST 


0.44 + 





.02 









,49 + 


,01 





,60 + 


,01 














L146 


IRAM 


5.28 ± 





,02 


5,28 + 


,02 


5 


,18 + 


01 





,39 + 


,01 





,25 + 





,06 





,25 + 0,06 




SEST 


5.22 + 





.02 






5 


,15 + 


,01 





,56 + 


,01 














L483 


IRAM 


4.95 ± 





,02 


5,13 + 


,02 


5 


,40 + 


,01 





,47+ 


,01 


-0 


,95 + 





,05 


-0 


,59 + 0,05 




SEST 


5.12 ± 





,03 






5 


,40 + 


,01 





,58+ 


,01 














Serp-S68N 


IRAM 


6.67 + 





,02 


7,83 + 


,02 


8 


,79 + 


,01 


1 


,32 + 


,04 


-1 


,61 + 





,05 


-0 


,73 + 0,03 


Serp-FIRSl 


IRAM 


9.64 ± 





,02 


7,92 + 


,04 


8 


,47+ 


,01 


1 


,11 + 


,01 


1 


,05 + 





,02 


-0 


,50 + 0,04 




SEST 


7.91 ± 





,04 






8 


,40 + 


,01 


1 


,29 + 


,01 














Serp-SMM5 


IRAM 


7.93 + 





,02 


7,84 + 


,02 


8 


,44 + 


,01 





,69 + 


,01 


-0 


,74 + 





,03 


-0 


,87 + 0,03 




SEST 


7.85 ± 





,02 






8 


,42 + 


,01 





,99 + 


,01 














Serp-SMM4 


IRAM 


7.35 ± 





,02 


7,69 + 


,02 


8 


,08 + 


,01 


1 


,84+ 


,03 


-0 


,40 + 





,01 


-0 


,21 + 0,01 




SEST 


7.59 + 





,02 






8 


,12 + 


,01 


1 


,80 + 


,02 














Serp-SMM3 


IRAM 


7.60 + 





,02 


7,80 + 


,02 


7 


,64 + 


,01 





,82 + 


,01 


-0 


,05 + 





,03 





,19 + 0,03 


18331-0035 


IRAM 


10.97 + 





,02 


10,99 + 


,02 


10 


,83 + 


,01 





,59 + 


,01 





,22 + 





,04 





,25 + 0,04 




SEST 


10.92 + 





,02 






10 


,84 + 


,01 





,71 + 


,02 














CrA-IRAS32 


SEST 


5.66 + 





,02 


5,58 + 


,02 


5 


,70 + 


,01 





,52 + 


,01 


-0 


,07 + 





,04 


-0 


,22 + 0,04 


L723 


IRAM 


11.25 + 





,02 


11,27+ 


,02 


11 


,14 + 


,01 





,63+ 


,02 





,17 + 





,04 





,21 + 0,04 


L673A 


IRAM 


6.94 + 





,02 


6,97 + 


,02 


6 


,92 + 


,01 





,47+ 


,01 





,03 + 





,05 





,11 + 0,05 


B335 


IRAM 


8.13 + 





,02 


8,14 + 


,02 


8 


,35 + 


,01 





,48 + 


,01 


-0 


,46 + 





,05 


-0 


,43 + 0,05 




SEST 


8.12 + 





,02 






8 


,35 + 


,01 





,53+ 


,01 














L1152 


IRAM 


2.60 + 





,02 


2,63 + 


,02 


2 


,70 + 


,01 





,40 + 


,01 


-0 


,24 + 





,06 


-0 


,19 + 0,06 


L1157 


IRAM 


2.46 + 





,02 


2,43 + 


,02 


2 


,69 + 


,01 





,71 + 


,01 


-0 


,33 + 





,03 


-0 


,36 + 0,03 


L1172 


IRAM 


2.73 + 





,04 






2 


,78 + 


,01 





,53+ 


,01 


-0 


,10 + 





,08 






L1251A 


IRAM 


-5.21 + 





,02 


-4,80 + 


,03 


-5 


,02 + 


,02 





,69 + 


,09 


-0 


,28 + 





,05 





,31 + 0,07 


L1251B 


IRAM 


-4.26 + 





,02 


-4,22 + 


,02 


-3 


,53 + 


,01 


1 


,00 + 


,01 


-0 


,73 + 





,02 


-0 


,69 + 0,02 


L1262 


IRAM 


4.49 + 





,02 


4,00 + 


,02 


4 


,21 + 


,01 





,45 + 


,01 





,62 + 





,05 


-0 


,46 + 0,05 



^The CS spectra were obtained at the Haystack Observatory with the exception of CrA-IRAS32, which was obtained at the SEST 
^The error we quote is the larger of one half the channel width and the formal la error from a gaussian fit to the line peak. 
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Table 4. H2CO 2o2 - loi and C^^S 2 - 1 Spectral Line Velocities 



Source 


Telescope 


VH2C0 


Vcsis 


AVc'345 






km 


km 


km 


L1527 


IRAM 




5.94 ± 0.02 


0.30 ±0.05 


13036-7644 


SEST 


3.42 ±0.02 


3.72 ±0.04 


0.53 ±0.10 


VLA1623 


SEST 


3.54 ±0.02 


3.58 ±0.05 


1.37 ±0.12 


WL22 


SEST 


3.58 ± 0.02 


3.87 ±0.02 


0.74 ± 0.04 


16293-2422 


SEST 


3.22 ±0.02 


4.02 ± 0.06 


1.78 ±0.21 


L483 


SEST 


5.31 ±0.05 


5.40 ±0.02 


0.71 ±0.05 


Serp-FIRSl 


SEST 


7.96 ±0.02 


8.53 ±0.03 


1.32 ± 0.08 


Serp-SMM5 


SEST 


7.87 ±0.07 


8.40 ± 0.03 


0.97 ±0.09 


Serp-SMM4 


SEST 


7.55 ±0.02 


8.11 ±0.05 


2.07 ±0.21 


18331-0035 


SEST 


10.82 ±0.02 


10.86 ±0.03 


0.37 ±0.07 


CrA IRS32 


SEST 


5.69 ±0.02 


5.72 ±0.05 


0.72 ±0.11 


B335 


SEST 


8.16 ±0.02 


8.43 ± 0.04 


0.59 ±0.10 


L1251B 


IRAM 




-3.74 ± 0.05 


1.42 ±0.12 



-27- 



Table 5. Statistical properties of the distributions of 6V 



Optically thick line 


Sample 


N 




No 


N+ 


mean ± s.e.m. 




E" 


H2CO 


all 


47 


15 


27 


5 


-0.14 ± 0.08 


0.06 


0.21 




Tboi < 70 K 


23 


12 


8 


3 


-0.28 ±0.13 


0.05 


0.39 




70 < Tboi < 200 K 


24 


3 


19 


2 


-0.02 ± 0.07 


0.81 


0.04 


CS 


all 


37 


14 


19 


4 


-0.14 ±0.07 


0.04 


0.27 




Tbol < 70 K 


19 


10 


9 





-0.28 ±0.10 


0.008 


0.53 




70 < Tbol < 200 K 


18 


4 


10 


4 


0.00 ± 0.09 


0.99 


0.00 



^N-, No, N+ arc the number of sources in the subsample with normalized velocity difference 5V < —0.25, 
-0.25 <SV < 0.25, and 0.25 < 6V respectively 

''Probability of drawing the sample from a zero-mean normal parent distribution, based on a Student's t-distribution 

"The blue excess is defined as E = (N- - N+)/N 
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Table 6. Line Asymmetries in Sources Having \SV\ > 0.25 in at Least One Line 



Source* 








Line Asymmetry'' 






K 


CS 2 - 1 


H2CO 2l2 


- Ill HCO+ 3-2 


HCO+ 4-3 


L1448mm 


56 


b 


r 


r 


r 


03256+3055 


74 


b 


n 






N1333-2 


52 


r 


r 




b(i) 


N1333-4A 


34 


b 


b 


b 


b 


JN1333-4B 


36 


b 


b 


b 


b 


L1551-5 


95 


n 


r 






L1551NE 


75 


r 


r 


b 




L1527 


60 


n 


b 


b 


b 


13036-7644 


63 


b(2) 


b 






VLA1623 


30 


b 


b 


n 


n 


WL22 


86 


b 


b 






16293-2422 


43 


b 


b 


b 


b 


L43 


79 


n 


r 






L483 


50 


b 


b 


r 


r 


S68N 


40 


n 


b 






FIRSl 


51 


b 


r 


r 


b 


SMM5 


59 


b 


b 






SMM4 


43 


b 


b 


b 


b 


B335 


37 


b 


b 


b 


b 


L1157 


62 


b 


b 


b 


n 


L1251A 


108 


b 


r 






L1251B 


91 


b 


b 






L1262 


104 


r 


b 







*The best cases for infall are indicated in bold face 

''The line asymmetry is designated b for SV < —0.25, n for —0.25 < 5V < 0.25, and r for 5V > 0.25. Except where 
noted, the CS and H2CO data come from this paper, and the HCO"*" line data come from Gregersen et al. (1997) 

References. — (1) Ward-Thompson et al 1996, (2) Lehtinen 1997 
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Fig. 1. — Observed spectral line profiles toward the sources with T;,,,; < 70 K in our sample ("Class 
0"). N2H+ 101 - 012 is drawn at the bottom of each box, H2CO 2i2 - In next, and CS 2 - 1 at 
the top whenever it was observed. The vertical scale is in main beam brightness temperature units. 
The horizontal scale has a range of 10 km s~^ in all sources. The dashed line indicates the N2H"'" 
velocity found with the hyperfine structure fit (the other six hyperfine components lie off the box). 
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Fig. 2. — Observed spectral line profiles toward the sources with T^oi > 70 K in our sample ("Class 
I"). The display format is the same as in Fig. 1. 
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Fig. 3. — Distribution of the line peak velocity difference 6V, defined in eq. [T|, for the full sample 
using: H2CO (left), and CS {right) as the optically thick line tracer. Note that more sources have 
negative (blue) than positive (red) values of 6V. 
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Fig. 4. — Relation between the velocity difference 5V based on H2CO (top)and CS (bottom) as 
the optically thick line, and bolometric temperature (T{,„;). Both plots show that most sources 
have blue-shifted values of 6V at low T^oj (Class sources). At higher T^oi (Class I sources) 
there are equal number of blue- and red-shifted sources. The uncertainty in 5V is typically ~0.05 
km s~^ (Table and the uncertainty in Tboi is typically of order 20 K (see text). 
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Fig. 5. — Relation between 5V determined using CS 2 — 1, and 5V determined using H2CO 2i2 — In 
as the optically thick lines; and using N2H"*" 101 — 012 in each case as the optically thin line. Note 
the good correlation in most sources. The uncertainty in 6V is typically ~0.05 km s^^ (Table |3[), 
comparable to the symbol size. 



